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Abstract 

The short distance QCD corrections to 6 — > 577 are calculated in the leading 

logarithmic approximation. The equivalence of operator basis reduction for 

S-matrix elements by using the equations of motion or by proving a low energy 

theorem is discussed. We apply the above results to the exclusive Bs 77 

decay. The branching ratio of this decay is found to be 5 x 10~^ in the 

Standard Model. We also found that QCD corrections modify considerably 

the ratio between CP-even and CP-odd two-photon amplitudes. 
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Rare B decays provide useful probes to the flavor structures of the Standard Model(SM) 
and its various extensions. Among them, the radiative decay B —>■ X^'y [|T| has recently 
received the most attention since it is not only within the reach of current experiments, 
where CLEO gives an updated branching ratio 

Br{B X,7) = (2.32 ± 0.51 ± 0.29 ± 0.32) x 10"^ (1) 

but also theoretically it has been calculated to the NLO accuracy i,i with the result 

Br{B X,7) = (3.28 ± 0.33) x 10-^ (2) 

in the SM. While this decay mode will be continually studied for some time, the upcoming 
B factories at KEK, SLAG and hadronic B projects at HERA, Tevatron and LHC offer 
new opportunities to study many more rare decay modes. With these facilities, it will be 
possible to study decay modes with branching ratios as small as 10^*^. In this work, we 
focus on the exclusive decay Bg '~f'~f, which, in spite of its small branching ratio, has a 
very clear signal where two monochromatic energetic photons are produced precisely back- 
to-back in the rest frame of Bg. Within the SM, the electroweak contributions to this 
process was calculated without QCD corrections some time ago The branching ratio 

so obtained was about 10~^. The calculation was later extended ||^ to the case of two-Higgs 
doublet models(2HDM). It was shown that Br{Bs 77) might be enhanced to 10~^ in 
certain regions of the parameter space. The current experimental bound on this decay is 
BriBs 77) < 1.48 X 10"^ 0. 

Before one goes on to study other new physics which potentially can influence this decay, 
it stands to reason to improve upon previous calculations by incorporating short distance 
QCD running from the scale mw to the scale m?,. This will be achieved by constructing 
an effective Hamiltonian for the quark level process b s'-ff at the scale fi = rrih. The 
amplitude for Bg 77 can then be determined once a method is devised for computing the 
hadronic matrix elements. 

At present, the effective Hamiltonian Heff{b 57) for 6 — 57 is already well understood 
T0| . It is natural to ask what is the corresponding effective Hamiltonian for b s'-ff. Is it 



identical to Heff{h 57) or more complicated with extra operators? The answer is that, 
to the order Heff{b — > 57) is sufficient to generate all contributions to b —* S77. 



This can be understood by either appplying the equations of motion ]T2[ or by applying an 
extension of Low's low energy theorem In the following we shall illustrate that the 

above two approaches are in fact equivalent. To proceed, we write the effective Hamiltonian 
for 6 ^ 57 as fill 



H^ffib .7) = -%v,:r,,^a(/i)o.(/i), 

i=l 



(3) 



where 
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07 = 


g 

-^Sia^^inisPL + mbPR)biFf,i, 


08 = 





(4) 

We remark that the decay b 577 also receives contributions induced by 7 and penguin 
diagrams, which are higher-order in electroweak couplings. Effects of such diagrams can be 
summarized into electroweak penguin operators which have similar structures as strong 



penguins O3 — Oq, and an 0(aem)-correction to the Wilson coefficient C3. It is well known 
that [|r^ penguin diagram contains the 0(m^ )-contribution which is not much suppressed 
compared to contributions of strong penguins. On the other hand, as we shall see later, the 
dominant contributions to 6 ^ S77 actually arise from Oi, O2 and O7. Therefore, in view 
of simplifying our analysis, we shall neglect electroweak penguin contributions. 

Concerning the effective Hamiltonian Hf.ff{b —>■ 57), we note that operator mixings 
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between the set, Oi, ■ ■ ■ ,0q, and the magnetic penguin O7, Og are generally regularization- 
scheme dependent. However, as pointed out by Ciuchini et al. in Ref. [1], the one-loop 
matrix elements of both b —>■ sj and b ^ sg arising from Oi, ■ ■ ■ , carry compensating 
scheme-dependence such that the total physical amplitudes are independent of regularization 
schemes. In view of this, one defines the so called "effective Wilson coefficients" |T^, C^^^ifi) 



and Cl"{^^), for operators Or and Og, for which the RG running is scheme independent. 
In 't Hooft-Veltman's (HV) scheme for 75, these effective coefficients coincide with original 
ones. Furthermore, in this scheme, only O7 contributes to the one-loop matrix element of 



6 — *• 57 |]T5[. Since we shall adopt the HV scheme for subsequent calculations, we drop the 
superscript "efF' on both Cj and Cg from now on. 

To see whether Oi, ■ ■ ■ , O7 are sufficient to generate all contributions to 6 — > 577, one 
notes that operators containing two photon fields are missing from Eq. (4). Photon fields 
which constitute operators in this category can come in 3 ways: from covariant derivatives 
D^, from field strength tensors F^^ or from some combinations of the two. It has been shown 
by Grinstein et al. in Ref. [1] that such operators can be eliminated or simply reduced to 
those in Eq. (4) by applying equations of motion. Therefore, H^ffib — » S7) is also the 
effective Hamiltonian for b 577. At this point, we wish to emphasize that the above 
simplifications by using the equations of motion are applicable only with respect to the S 
matrix elements ||12|. The Green functions depend on the eliminated operators, however. 

The above prodecure based on applying the equations of motion has a close link with 
the low energy theorem derived by us some time ago ||^. Let us illustrate this with the 
operator = sip ^bi [|1^]. Upon applying the equation of motion for the b quark, it is 
easily seen that Q4 = imfsLbR. In this new form, Q4 gives rise to the flavor changing self 
energy. However this contribution will be subtracted out by the on-shell renormalization 



sS(p2 = ml) = 0, S(P2 = rnl)b = 0. (5) 
Hence, with equations of motion one concludes that Q4 gives no contributions to 6 — > S77. 



This fact is also realized by the low energy theorem. There one takes Q4 as it is initially 
defined and computes Q^s contributions to b — s self energy, 3-point b S'y and 4-point 
b —>■ S77 vertices. Its total contributions to 6 ^ 577 are given by summing up the irreducible 
and reducible diagrams where the latter come from attaching bremsstrahlung photons to 
the fiavor changing self-energy and the three-point b ^ S'y vertex. The renormalization 
condition, Eq. (5), ensures that the end result will be local, and one can further show that 
the extraneous contributions cancel out due to gauge invariance p. Such a cancellation 
among reducible and irreducible digrams is what was referred to as the generalization of 
Low's low energy theorem It is clear that this theorem is essentially a verification of 

general arguments based upon applying the equations of motion. 

With the above issues clarified, we are ready to compute b — > S77. The amplitude can 
be separated into irreducible and reducible parts, as shown in Figs. 1 and 2. 

Figure 1 



Figure 2 

We begin with irreducible diagrams which are more involved. There is only one type of 
diagrams which are depicted generically by Figure 1, where the four-point vertex is a result 
of inserting 4-quark operators Oi — Oq. While in the case of Oi and O2, Figure 1 consists 
of only one diagram, there are, however, two different ways of contracting in the case of 
O3 — Oq, as the light quark field q appearing in these operators can be external particles as 
well. To simplify the algebra, we apply Fierz rearrangements to four-quark operators when 
necessary. This rearrangement is legitimate because the 75 we use is simply the product of 
7°, ■ ■ ■ , 7^. The Fierz identities we have used are 
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{qiq2)v-AiQ3q4)v-A = {qiqi)v-A{'im)v-A^ 
{qiq2)v-A {q-m)v+A = -2 + 75)^4) (53(1 - 75)^2) • 

The amplitude for the irreducible diagrams can be written as 
Mir = Vt,Vtb X 



where 



971 



u{p' 



J2 AgJ{ml)YPLR^.up + (msK{ml)PL + m,K{ml)Pa) T, 



L q 



+C {-msL{ml)PL + m,L{ml)Pn) e^^^pk'^k^] u{p)e''{h)e''{k2 



(6) 



(7) 



P 



+ ki ■ k2e^,yp„{k2 - kiY 



(8) 



is a third rank tensor originally constructed by Rosenberg and Adler [18|; and 



T^u — (^2,^^1,1/ — g^uki ■ ^2) 



(9) 



Here we adopt the convention: £0123 = 1- The coefficients A^'s, B and C are linear combi- 
nations of Wilson coefficients and are given by 

= (C3 - C^)N, + (C4 - Ce) 

Ad = \{{C^-C,)N, + {C^-C^)) 

A, = ((Ci + C3 - C^)N, + (C2 + C4 - Ce)) 
1 



A, = Ah 



((C3 + C,- C,)N, + (C3 + C,- Ce)) 



B = C = --iCeN, + C,). 



(10) 



Note that the above Wilson coefficients are evaluated at scale /x = rrih. Finally the fuctions 
J{m?), K{m?) and L{m^) are defined by 

J(m2) = Jn(m2), 
ir(m2) = (4/n(m2)-/oo(m2)), 



(11) 



with 

Ipg{m') = I dx I dy „ , -. (12) 



1 rl-x x'Py'^ 



JO 711? — 2xyki ■ k2 — ie 
For the reducible diagrams, the one loop contributions with Oi-Oq insertion actually 
vanish, similar to what occurs in 6 ^ 57 in HV-scheme. Therefore the amplitude for 
reducible diagrams is identical to its lowest order form p except for replacing C7{'mw) by 
Ci^rrif,). Hence 

OTT 



^-^)(r,pa,^k^,k^ + 2i{^ \)o.pk^2 



\_ p ■ ki p ■ K2 p ■ ki p ■ ki 

■{nn,PR + m,PL)u{p)e^'{ki)e''{k2) + k^, fi< — >k2,iy. (13) 

From Eqs. (7) and (13), we obtain the total amplitude for b —>■ sj'-f. 

To calculate Bg — > 77, one may follow a perturbative QCD approach which includes a 
proof of factorization, showing that soft gluon effects can be factorized into Bg meson wave 
function; and a systematic way of resumming large logarithms due to "hard gluons" with 
energies between 1 GeV and mf,. Since such an approach is technically rather demanding with 
uncertain reliability, we shall adapt a phenomenological approach where the long distance 
effects are replaced by a few non-perturbative parameters. In other words we simply evaluate 
the hadronic matrix element of Mipi + Mji, relying on a phenomenological model. Before 
doing so, it is important to note that Mji is apparently non-local due to internal 6 or s quark 
propagators. To handle these non-local terms, one observes that the b quark inside the Bg 
meson carries most of the meson energy, and its four velocity can be treated as equal to that 
of Bg. Hence one may write b quark momentum as p = rribV where v is the common four 
velocity of b and Bg. With this parametrization, we have 

p- ki = rubv ■ ki = IrUbmB, = p ■ k2, 
p' ■ki = {p-ki- k2) ■ ki = -^uibX^Bs - ^b) = p' ■ ^2, (14) 

where the second equation is based on a constituent picture that b and s quarks share the 



total energy of Bj^. Therefore —p' should be taken as the four momentum of a constituent 
s quark. With Eq. (14), is readily made local. We then compute the amplitude for 
Bs 77 using the following relations 

(0|S7^75&|5.(P)) = -^fB.P,, 

{Q\sj,bmP)) = tfs^MB, (15) 

where /^^ is the Bg meson decay constant which is about 200 MeV according to recent 
Lattice QCD calculations ]TP[ ]. 

The total amplitude is now separated into a CP-even and a CP-odd part 

T{B, ^ 77) = M+F^,F^- + tM-F^,F>^r (16) 

We find that 



and 



= - ' fB.^ByLyt^> \Bm,K{mi) + ^ ) , (17) 



M- = ^^^^fBmsy:,V,, (e^b A,J(mJ) + m,BL{ml) + , (18) 



where A = rriB, — mt,. The decay width for Bg '~f'~f is simply 



T{Bg 77) = J^i\M^\ + |M+r). (19) 



To obtain numerical results, we have set light quark masses to zero and used [^] rrit = 
175 GeV, rrif, = 4.8 GeV and rric = 1.5 GeV. Furthermore, we take = 5.37 GeV, 
a = and Vj*Vtb = 4 x 10~^. The numerical values for Wilson coefficients Ci — Cg 
evaluated at /i = rrib are listed in Table I. With the above input parameters, we find T{Bs — *■ 
77) = 2.0 X 10~^° eV which amounts to a branching ratio Br{Bs — * 77) = 5.0 x 10^'^, 
for the given F^*"' = 4 x 10~^ eV. If QCD corrections are not included, namely taking 



^Note that the momentum of s quark is —p' as p' denotes the momentum of s quark in 6 ^ 577. 
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C2 = C2{mw) = 1, C7 = CT^mw) and setting all the other Wilson coefficients to zero, one 
obtains T{Bs — > 77) = 1.3 x 10^^'^ eV. Therefore including QCD effects has increased the 
rate of Bs —>■ '~f'~f by more than 50%. 

It is interesting to note that QCD correction modifies contributions of irreducible and 
reducible diagrams in opposite ways. It enhances the contributions of reducble diagrams 
through the enhancement of C7, namely C7{mb)/C'j{mw) ~ 1-5. Such an enhancement 
is already well known in the decay b 57. As a contradistinction, QCD suppresses 
the contributions of irreducible diagrams since its effect essentially replaces C2{mw) with 
NcCi{mi,) + C2{mi))^. The latter is much smaller than the former due to cancellations be- 
tween NcCi{mh) and C2(mfe). Because QCD gives distinct effects to Mjr and Mr, the 
relative magnitude of M"*" and M~ is also modified accordingly and drastically by QCD. 
From Eqs. (10), (17) , (18) and Table 1, it is clear that the magnitudes of and M~ 
are almost identical, both dominated by Cj but somewhat corrected by the suppressed 
NcCi{mb) + C2(mfc). Numerically we have |M"*"|^/|M^|^ = 0.80. The corresponding ratio 
without QCD corrections is 0.38 which is twice smaller. 

On the experimental side, as mentioned before, the future B factories are capable of 
probing B decays with branching ratio as low as 10^^. Therefore one expects the decay 
Bs 77 to be seen in these future facilities. It is, however, more challenging to separate 
partial amplitudes M"*" and M~ because it will require measuring the spin correlation of 
final state photons. 

At this point, we wish to state the obvious interdependence between 77 and 

inclusive B Xs'j decays. It is well known that the later decay depends on Cjlmw) at the 
tree level, and Cs{mw) and C2{mw) through QCD-induced operator mixings ||l|. CLEO 



measures some combinations of these coefficients, and thus imposes constraints |1^J21 



^Since C3{mh), • • • , C^{mh) are much smaller than d(mf,) and C2{mb), the dominant contribution 
to Mjfi is proportional to Ac which is approximately equal to NcCi{mf,) + C2{mi,). 
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on various extensions of the Standard Model, such as 2HDM or Minimal Supersymmetric 
Standard Model(MSSM). Since Bg 77 depends on the same set of Wilson coefficients, its 
sensitivity to physics beyond the Standard Model complements the corresponding sensitivity 

Before closing, we like to comment on the inclusive B ^s77 decay. This process has 
recently been explored without the inclusion of QCD corrections |2^. The branching ratio 
is found to be 10"'' in the SM. Experimentally, due to the appearance of hadronic final state, 
it seems to be more difficult to analyze this process than one can do for Bs —* 77. Besides, 
photons in this decay may be soft or collinear which may require extra treatment such as 
the necessity of imposing kinematical cuts. On the theoretical side, there exists a further 
complication in that photons can be emitted from the spectator quark which forms B with 
b. Such spectator effects, while neglected in Ref. [23], need to be studied carefully in order 
that additiona 

1 informations on C2 and Cj can be extracted from B —>■ ^^77. 

In conclusion, we have outlined the equivalence of procedure by using the equations of 
motion [|T^] and that by the low energy theorem P,p!B[ to reduce the operator basis, as they 
are applied to 6 — 577. We have also computed the exclusive decay Bg — > 77 with the 
effective Hamiltonian Heff{b 57). The branching ratio is found to be 5 x 10^^ which 
is 50% larger than the result without QCD corrections. Finally we have argued that one 
should include spectator effects in the study of inclusive B — > ^^77 decay. 

Note added. After completing this work, we became aware of a paper by G. Hiller and E. 
O. Iltan, [hep-ph/9704385| , where the QCD correction to Bg '~f'~f is also discussed. However, 
the crucial contribution from Oi and other contributions from O3, ■ ■ ■ ,0q are not taken into 
account in that paper. 
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FIGURES 

FIG. 1. The generic irreducible diagram contributing to b ^ 577. The cross in a circle denotes 
insertions of operators, Oi, • • • , Oq. We omit the diagram with two photon lines interchanged. 

FIG. 2. The reducible diagrams contributing to b ^ S77 where the interchange of two photon 
lines is assumed. The cross in a circle denotes the insertion of operator O7. 
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TABLES 

TABLE I. Wilson coefficients Ci{iJ,) at fi = mi, = 4.8 GeV in tfie leading logarithmic approxi- 
mations. These values are obtained by taking Mw = 80.2 GeV and as{Mz) = 0.117 



Ci(m) 






C4(/i) 


C5(/^) 




C7(m) 


C8(/X) 


-0.222 


1.09 


0.010 


-0.023 


0.007 


-0.028 


-0.301 


-0.144 



15 




Figure 2 



